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ABSTRACT 

Using Herschel's HIFI instrument we have observed [Cn] along a cut through S140 and high-y transitions of CO and HCO^ at two positions 
on the cut, corresponding to the externally irradiated ionization front and the embedded massive star forming core IRSl. The HIFI data were 
combined with available ground-based observations and modeled using the KOSMA-r model for photon dominated regions. Here we derive the 
physical conditions in S140 and in particular the origin of [Cii] emission around IRSl. We identify three distinct regions of [Cii] emission from 
the cut, one close to the embedded source IRSl, one associated with the ionization front and one further into the cloud. The line emission can 
be understood in terms of a clumpy model of photon-dominated regions. At the position of IRSl, we identify at least two distinct components 
contributing to the [C ii] emission, one of them a small, hot component, which can possibly be identified with the irradiated outflow walls. This is 
consistent with the fact that the [C ii] peak at IRSl coincides with shocked H2 emission at the edges of the outflow cavity. We note that previously 
available observations of IRSl can be well reproduced by a single-component KOSMA-t model. Thus it is HIFI's unprecedented spatial and 
spectral resolution, as well as its sensitivity which has allowed us to uncover an additional hot gas component in the S140 region. 

Key words. ISM: structure ~ ISM: kinematics and dynamics ~ ISM: molecules -PDRs - Submillimeter 



' 1. Introduction 

Massive star formation poses a challenge to the observer, since 
. massive stars are formed predominantly in clustered environ- 

■ ments and exhibit a number of entangled energetic phenomena in 
' a small region. The S140 region shows several phenomena asso- 
ciated with (massive) star formation, such as outflows and strong 
UV irradiation both from internal and external heating sources 

, creating photon dominated regions (PDRs). Thus, S140 could 

■ also act as a template for extragalactic star formation, which is 
' traced by strong fine structure lines stemming from PDRs. 

\ At a distance of 910 pc, the BOV star HD211880 ion- 
' izes the edge of the molecular cloud L1204, cr e ating SI 40, 
] a visible Hii region and PDR dCrampton et al.L Il974l) . The 
dense molecular cloud hosts a cluster of embedded mas- 
sive young stellar obj ects (YSOs) only 75" fro n t the H 11 re- 
gion (e.g. Beichman etlil [T979HMinchin et"aLL Il993h . S140 
and its surroun dings were extensively ob s erved in molecu- 
lar tracers (e.g. 'Minchin et all I1993L 1 1995b Ivan der Tak et all 
[2OOO; Persson et al., 2009) and the fine structure line s , [Cn ] 
(Boreikoetal., 1990; White & Padman, 1991; Li et al., 2002), 
[Ci] (.Minchin etaL .1994) and [Oi] (Emerv et aL 1996.). An 
ionize d co mponent aroun d IRSl has been detected by iHoard 
(I2OO6') and 'Trinidad et al.' ('2007') at high resolution in the cm 
continuum observations. Spaans & van Dishoeck ( 1997) model 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



SI 40 with a nearly edge on geometry, a clumpy, inhomogeneous 
medium and a radiation field of 150 Drainq^ units. 

Two bi-polar outflows can be found in t he S140 cloud 
(e.g. 'Wei gelt et all l2002t iPreibisch et"dl l2002h . One is traced 
through observations of bow-shock like features and H2 knots, 
indicating hot, shocked gas, lined up northeast of IRS 1 over sev- 
eral degrees, while the other, perpendicular to it, has been ob- 
served in CO. 

In this paper, we describe and model the [Cii], high-7 CO 
and HCO^ observations in the S140 region, observed as part 
of the Warm and dense ISM key proje ct (WADI, Ossenkopf e t 
al., KPGT_vossenko_l) using the HIFI (Ide Graa uw et al.l |2010h 
instru ment on the Herschel Space Observatory (Pilbrat t et all 
I2OIO V The goal of these studies it to investigate and compare 
the physical properties of the different gas phases in two regions 
of the S140 cloud, the star forming region IRSl and the ioniza- 
tion front (IF from now on). 

2. Observations & Analysis 

2.1. HIFI observations 

For this paper, two single pointing, frequency switching spec- 
tra in bands la, 4a and 4b from the Priority Science Phase 
and a [Cii] OTF map from the Science Demonstration Phase 
were analyzed. The single pointing data were taken towards 
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S140-IRS1 at RA=22hl9ml8.21s, Dec=63°18'46.9" (J2000) 
and towards a position in the IF at RA=22hl9ml 1.53s, 
Dec=63°17'46.9" (J2000). The fully sampled [Cii] map is cen- 
tered on RA=22hl9ml8.88s, Dec=63°18'52.9" (J2000), has a 
length of 206" and crosses the PDR region from the IF through 
IRS 1 into the dense molecular cloud (Fig.[Tl right). The data pre- 
sented here were taken with the wide band spectrometer (WBS) 
at 1 . 1 MHz resolution. For the FS W data, the upper and lower 
sidebands were processed separately. All observational param- 
eters are summarized in T able [1] The data reduction was done 
using HIPE 3.0 (IOttetaUl2010l) . A forward efficiency of 96% 
and main beam efficiencies of 0.72, 0.74, 0.76 and 0.65 at 500, 
950, 1112, 1900 GHz, respectively (R. Moreno, priv. comm.), 
were used to convert ffie data to main beam temperatures r^b- 

2.2. Complementary data 

We have combined the HIFI data with observations of all CO 
isotopes in the 2-1 and 3-2 transitions, of '^CO in the 4-3 and 
7-6 transitions, HCO+ and H'^CO^ in the 3-2 and 4-3 transi- 
tions, and both fine structure lines of atomic carbon [CI], taken 
at the KOSMA 3m sub-mm telescope, the IRAM 30 m tele- 
scope, and the JCMT. For easy comparison of flux values, all 
data have been smoothed to the 80" resolution of the KOSMA 
345 GHZ data. The procedure to derive appropriate scaling fac- 
tors for the single-point observations with different beams, based 
on the spatial distribution measured in the SPIRE ISOj-im and 
SCUBA 450//m map, is described in detail in the Appendix. 

2.3. PDR modeling 

To model ffie physical parameters in IRSl and the IF, w e fit 
the data by the KOSMA-t PDR code (iRollig et al.Ll2006l) . The 
PDR is represented by an ensemb le of spherica l clump s with a 
power-law clump-mass spectrum dCubick et all l2008l) . Details 
are summarized in Appendix B. Each clumpy PDR ensemble 
has five free parameters: the average density of the clumps 
in the ensemble, n^ns, the ensemble mass, Mp„5, the UV field 
strength, x given in units of the Draine field, and the minimum 
and maximum mass of the clump ensemble, [M„„„,M„,a,:]. We 
fit absolute line intensities, using the available ground-based 
observations and the HIFI lines of the CO isotopes, HCO^, 
atomic and ionized carbon. Simulated Annealing (Ingber, 
11993:) was used to find the optimum parameter combination. 
Our model ignores mutual shielding and illumination effects 
between different clumps and/or ensembles, i.e. optical depth 
effects are only considered within individual clumps . 



3. Results 

3.1. [Cii] cut 

Fig. [1] left, shows the measured position-velocity diagram of the 
[C n] cut. We can distinguish three peaks. The first peak (IF-peak 
from now on), at around -1.2', is located right at the edge of the 
rim of the IF as seen by White & Padman (1991) and represents 
the [C n] emission from the IF. A second, much broader peak is 
found around 0', in the vicinity of IRSl (therefore called IRSl- 
peak). The third peak is located at 1' into the dense cloud. 

At the IF position, we see symmetric [C ii] line profiles and 
no gradient in line width or centroid velocity. We find neither in- 
dications of a pressure gradient nor of an evaporation flow from 
the PDR surface. Looking at archival IRAC 8 //m data, the peak 
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Fig. 1. Fig. 1 Left: Position-velocity diagram of the [C n] emis- 
sion from the IF into the cloud, centered on RA=22hl9ml8.88s, 
Dec=63°18'52.9" (J2000). IRSl is at found at -7.5". Right: 
IRAC 8jum map of S 140. Shown are the two positions in which 
individual HIFI measurements are taken as well as the location 
of the [C ii] cut. HD21 1880 is located to the left of the pv dia- 
gram, and to the bottom right of the cut in the right panel. 
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Fig. 2. HIFI spectra at IRS 1 and the IF. The additional ghost line 
at 4 km s"' in the HCO+(6-5) and '^CO(5^) data is an artifact 
from the FSW deconvolution of the two lines, close in IF fre- 
quency but stemming from different receiver sidebands. Some 
lines are scaled for a better comparison. The two dashed lines 
represent the Vki at -7.1 km s"' (IRSl) and -8 km s"' (IF). 



of the [C ii] integrated intensity at the IF-peak is consistent with 
that of IRAC 8 fim, tracing the PAH and evaporating VSG emis- 
sion. 

The emission at the IRSl -peak is 10" off'set from the young 
stellar object IRSl and coincident with the locations of the 
dee ply embedded source Submm-2, which was first described 
by Minchin et al.l d 1 9951) and the strong H2 knot #1 found 
by P reibisch et al.l (|2002|) . Compared to the K-band data of 



iWeigelt et al.l (120021) . the [C 11] peak is close to the edge of what 
they describe as the wall of the outflow cavity carved into the 
material by the jet. The strong [C 11] emission seen here is prob- 
ably created at this interface, either by irradiation of the out- 
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Table 1. Summary of the HIFI observations used and basic results 



Transition 


Vlinc 


HPBW 


Mode' 


S140-IRS1 


S140-1F 




Source 


nns"^ 


J r^''' 


I'Lsr 


Av 


^source 


rnis^ 


J ™^' 


I'lsr 


Al' 




[GHz] 


" 




(Si 


[Kl 


(Kkms-') 


(kms-') 


(kms"') 


[s] 


[K] 


(Kkms-') 


(kms"') 


(kms-') 


CII 


1900.537 


12 


OTF 


12 


0.9 


59.63 (0.99) 


-7.26 (0.05) 


6.82(0.15) 


12 


0.9 


96.91 (0.60) 


-8.00 (9.2e-03) 


3.05 (0.02) 


"C0(5^) 


550.926 


43 


FSW 


480 


0.02 


85.74 (2.48) 


-7.24 (0.06) 


4.02(0.14) 


240 


0.03 


4.56 (0.04) 


-7.89(0.01) 


2.85 (0.03) 


C"*0(5^) 


548.831 


43 


FSW 


480 


0.02 


17.29 (0.07) 


-7.24 (5.7e-03) 


2.98 (0.01) 


240 


0.03 


0.47 (0.03) 


-7.66 (0.08) 


2.21 (0.17) 


HCO+(6-5) 


535.062 


43 


FSW 


480 


0.02 


38.96 (4.90) 


-6.77 (0.23) 


3.75 (0.57) 


240 


0.16 


0.67 (0.28) 


-7.52 (0.74) 


3.48(1.63) 


C"*0(9-8) 


987.56 


23 


FSW 


66 


0.1 


2.60 (0.07) 


-7.03 (0.04) 


2.92 (0.09) 












"CO(10-9) 


1101.35 


20 


FSW 


156 


0.15 


23.41 (0.73) 


-6.99 (0.06) 


3.93(0.15) 













' OTF = On-The-Fly, FSW = frequency-switch, OFF position = 22hl8m37.01s, 63° 14' 17.9" 



' At a velocity resolution of Av=0.7 kms 



flow walls (e.g. iBruderer et al.L l2009allbh or shock interactions, 
as traced by the H2 knots. The [C n] lines are very broad with 
a pronounced red wing tracing an outflow of ionized carbon, 
which is visible in Fig. |2] between -5 km s"' and km s"'. A 
2-component Gaussian profile with components at -8 km s ' 
and -5.5 km s ' fits the observed [Cii] line profile. Since the 
-5.5 km s ' component of the [C n] emission is not observed in 
any of the other molecular tracers, it does not arise in the out- 
flows, but may originate in the ionized medium. Between 
the two [C 11] peaks, the [C 11] emission drops significantly, in- 
dicating a spatial separation between the IF component and the 
component around IRS 1 . This drop in [C n] intensity between 
the IF-peak and IRS 1 -peak roughly coincides with a local mini- 
mum of the dust temperature, where neither the heating from the 
internal IRS 1 source nor the external star (Lester et al., 1986) are 
efficient. 

In the region where the third peak is seen in [C 11], the SPIRE 
250 i-im map also shows a secondary peak (Juvela et al., priv. 
comm.), which points to a deeply embedded source. The [Cii] 
line has pronounced wings, at a slightly higher velocity than in 
IRS 1 . This might indicate a separate outflow or material shocked 
by the outflow from IRS 1 . 

Several spectra were taken with HIFI at the IRS 1 and IF po- 
sition. Overlay of those spectra of different tracers shows that 
with the exception of [C 11] all the molecular lines at IRS 1 have 
similar line shapes (Fig. |2]l. The [C 11] line shows a strong red 
wing, which was also not de tected in the p r evious coarser spatial 
resolution observations by IBoreiko et al.l (II990 '). The centroid 
velocity, Vcen, for all the molecular lines range from -6.8 to -7.4 
km s ' and agree well with a source veloci ty of -7.1 km s~' 
as derived from previous observations (e.g. ' van der Tak et al.L 
l2000t [Persson et al],l2009h . Although somewhat fainter in inten- 
sities (particularly HCO^ 6-5), several of the molecular lines are 
still detected at the IF position. The line widths at the IF position 
are up to a factor of two narrower compared to the widths at the 
position of IRS 1 . At the IF position the centroid velocities fall 
between -7.5 to -8.0 km s"' , corresponding to a slight blue shift 
compared to the velocities at the IRS 1 position. 
Detailed discussion of the [C 11] line profiles will be taken up in 
a later paper by Mookerjea et al. (in preparation). 

3.2. PDR model in IRS1 

Since the PDR region arising from the irradiated outflow walls 
is illuminated both from the outside by the BOV star HD211880 
and from the inside by the cluster of YSOs around IRS I, we 
assume a two component fit with a hot component originating 
at IRSl with strong FUV illumination, but only a small frac- 
tion of the total mass, and a cooler component that provides 
the bulk of the material. Owing to the much coarser resolution 
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Fig. 3. Clumpy PDR model fit to the observed CO and HCO+ 
line integrated intensities in IRS 1 , shown as function of the up- 
per level energy. HIFI measurements are depicted as open trian- 
gles, complementary data points as open circles. The range of 
model intensities spanned by the total mass range of 54-250 M© 
for the cool component is shown as shaded area for the 7-lines. 

Table 2. The parameters for the two ensembles in the PDR 
model in IRSl. 



<ne„s> 


Men. 


X 


[M,„i„,M„„] 


N" 


./v 


ri' 


(cm--') 


(Mo) 


(Draine) 


(Mo) 








Hot: 














1.8 10' 


14 


2.3 10' 


[0.008,7] 


151 


0.42-0.73 


1.3 


Cool: 














1.3 10' 


54-250 


27 


[0.008,27] 


413- 1913 


0.07 - 0.34 


0.6-2.6 



" Number of all clumps; ^' assuming hot clumps are situated approxi- 
mately 7" away from IRS 1 ; /v volume filling factor within a shell with 
^v,innci=5" and i?v,oiitci = l 1-13"; /a area filling factor withi?A=7" for the 
hot component. For the cold component, i?v.inner= 11-13", i?v,oiitcr=40" 
and R/, = 40" 



of the available hydrogen and carbon radio r ecombination lines 
(ISmirnov et al.LII995l;IWvrowski et"aniI997h it is diflicuh to as- 
certain the fraction (if any) of [C 11] emission being contributed 
by the ionized gas associated with IRS 1 . Based on the Gaussian 
decomposition described in Sect. 3.1, the [Cii] intensity from 
the PDR must fall between the 60 K km s ' of the total inte- 
grated line and the 37 K km s • of the -8 km s ' component 
only. Consequently, we obtain two different limiting fits that dif- 
fer mainly by the mass of the cool component falling between 
54 and 250 Mq. The best fit result is shown in Fig. [3] The cor- 
responding model parameters for both components can be found 
in Table |2] Figure [3] also indicates the contributions of the two 
components to the emission. 

The reduced x' of the fit is 2.4. We applied a generic 30% 
error to all data points. The model slightly underestimates the 
mid- J lines of the CO isotopologues and of HCO^. This is the 
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Table 3. PDR model fit to the observed atomic fine structure fine 
integrated intensities in IRS 1 . Shown are the ranges spanned by 
the two fits. 

species total hot cool obs. 

(Kkms-') (Kkms-') (Kkms-') (Kkms"') 



[Cii] 
[Ci](l-0) 

[Ci](2-1) 



19.3-21.9 
9.7-11.0 
6.9-8.5 



16.5 
1.5 
1.7 



2.8-5.4 
8.2-9.5 
5.2-6.8 



37. 

11.4 

20.5 



energetic region where neither the cool nor the hot component 
are very bright, so that this may hint towards a third intermedi- 
ate component. However, the available data do not allow to reli- 
ably derive parameters of a third component. To see the modeled 
FUV flux of the hot component, the clumps have to be at a dis- 
tance of 7" from IRS 1 . At this distance, an area filling factor on a 
spherical shell of 1.3 is obtained, i.e. the hot component shields 
the cool ensemble from most of the FUV flux. To obtain rea- 
sonable values of volume filling factors, i.e. 0.42 < /v < 0.73, 
the inner and outer radii of the shell have to be roughly equiva- 
lent to 5" and 1 1-13", respectively. While at the larger radii, the 
FUV flux seen by the hot component would be somewhat lower 
than predicted by our model, this discrepancy is still small com- 
pared to the uncertainties in the line estimates and thus the model 
parameters. The model parameters of the cool ensemble are par- 
ticularly interesting. The FUV field is lower than expected from 
the externally illuminating star HD211880. This suggests that 
the prominent edge-on PDR at the IF blocks the FUV radiation. 

3.3. PDR model in the IF 

For the IF position we are not yet able to provide a full model 
fit to the data because of the lower numbers of detected lines. 
We can, however, derive some qualitative conclusions. The mid- 
and high-y lines of CO and HCO^ can only be explained by 
a hot gas component with densities of the order of 10'' cm"^ 
and a local FUV intensity in the order of 10^ Draine units. The 
external radiation from HD 2 11 880 cannot provide more than a 
few thousand Draine units. An additional, so far unknown UV 
source seems to be present. The total gas mass in the beam at the 
IF is on the order of 10 Mq. 



4. Discussion and Conclusions 

About 90% of the [Cii] emission around IRSl stems from 
the hot component near the embedded IR sources. A different 
heating mechanism accounts for the [CII] from the emission in 
the ionizing front. This is reflected in our choice of two separate 
models for the two regions. 

Previous models of ^PDR emission in SI 40 

(ISpaans & van DishoeckL Il997h describe the S140 region 
as an inhomogeneous medium with (riu) = 7 x lO^cm"^, a filling 
factor of 30% and clumps of 0.04 pc size. The radiation field, 
X, of 150 Draine units they derive is much lower than the one 
derived in this work, however, these authors include the IRSl 
cluster as an embedded heating source in their model, yet do 
not distinguish between the two regions of [C ii] emission - the 
IRS 1 region and the IF. Moreover, they could not include high-7 
CO lines, which are important tracers of the hot gas. In our two 
component approach, the high-/ CO lines on one hand tell us 
how strong the hot ensemble has to be illuminated by the FUV 
and they determine how the total mass is distributed between 



hot and cool component. High-7 CO lines will be especially 
valuable since they, in the context of PDR models, can only be 
produced by the very high FUV fields while lower-7 lines have 
always more mixed origins and hence pose weaker constraints 
for the multi component model. With the exception of the 
high-7 '^CO(10-9) and C'**0(9-8) fines most other data points 
can be explained with much more moderate FUV fields. 

When repeating our calculations using only pre-HIFI data 
we reproduce the parameters obtained in prev ious models 
(e.g. ISpaans & van Dishoecld 119971: iKosteii Il998h . It becomes 
apparent that only HIFI data can reveal additional components 
deep inside the cloud. The low mass of the hot component 
together with the low FUV field of the cool component suggests 
that only a small fraction of the material sees the full FUV field, 
supposedly where the outflows created cavity walls at the edge 
of the ambient material. The bulk of the material is effectively 
shielded from the FUV photons. 

At the IF, first qualitative results hint at a so far unknown 
heating source and a two component density structure. A full 
PDR-model on this position (Rollig et al., in preparation) will 
reveal further properties of the material at the IF. 

HIFI's high spectral resolution observations allow us to 
identify the general velocity fields of the molecular and ionized 
gas around IRSl and trace the dynamics of the region, such as 
the blue outflow emission associated with the SPIRE-peak. This 
is possibly due to another, still deeply embedded object, which 
is driving its own outflow. 
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Appendix A: Comparison of line profiles taken with 
different beams 

The different lines discussed fiere were measured at various tele- 
scopes and at different frequencies so that they all represent 
somewhat different spatial resolutions. To allow a comparison in 
terms of a physical interpretation, they have to be translated to a 
common resolution, so that they stem from the same area on the 
sky. As a reference, we use 80", the resolution of the KOSMA 
observations in the 3-2 transition of the CO isotopes. In princi- 
ple, all data taken at a finer resolution can be resampled to this 
beam if the mapped area is large enough. 

Unfortunately, most HIFI observations were only single- 
point observations, not full maps, so that such a convolution is 
impossible. To derive scaling factors that describe the translation 
between the measured intensities and the intensity that would be 
obtained in an 80" beam, we have to assume a source geome- 
try of the emission. Instead of using any analytic geometry, we 
use the actually measured distribution of warm dust seen in the 
sub-mm. By assuming that the spatial distribution of all PDR 
tracers roughly follows the warm dust we derive scaling factors 
for the line intensities at different beam widths, by convolving 
the sub-mm continuum map with the different beam sizes and 
picking ratios between the convolved intensities at the measured 
positions. 

As a continuum map at IRS 1 , we used the combination of 
the SPIRE 250jum (18" beam siz^) and SCUBA 450//m im- 
age (8" beam size, Holland et al. 1999), because IRSl is satu- 
rated in SPIRE 250 //m, whereas the observed area of SCUBA 
450 ;um is too small to make a convolution map with the beam 
size of 80". We regrid the SCUBA 450/im image to the same 
grid as SPIRE 250 ;um, determine the scaling factor between the 
SCUBA 450yum and the SPIRE 250 yum maps from the overlap- 
ping area, and replace the saturated pixels of SPIRE 250;um with 
the scaled SCUBA data. As this combination implies some arbi- 
trariness, we have tested four different approaches. We either de- 
rive the scaling factor in a least-square fit using (1) all the valid 
overlapping pixels or (2) only overlapping pixels with SPIRE 
250;um > 100 Jy/beam. When replacing SPIRE pixels, we either 
replace only saturated pixels (1), or replace the full square area 
(5x5 pixels, 2) containing all saturated pixels. The combina- 
tion of these provide 4 different beam scaling factors, which are 
consistent within 3%. Taking the average of these 4 values, we 
derive the final factors as shown in Table lA.ll A direct convo- 
lution to 80" was possible for the ground-based maps that were 
observed with a smaller beam, such as the JCMT CO(3-2) map. 
All the resulting intensities are summarized in Table lA.2l 



Table A.l. The scaling factors between the different beam size 
observations, which should be multiplied to the line intensity to 
estimate the one in a 80" beam. 
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Fig. A.l. a) The spatial distribution of warm dust obtained by 
combining the SPIRE 250yum map with the SCUBA 450yum map 
(see text). The HIFI beam is assumed to "see" convolved images 
at beam sizes of (b) 42", (c) 55", and (d) 80". The color scale 
is common for all the figures. The coordinates are relative to the 
position of IRS 1 . Thus the flux ratios at (0,0) in these images are 
the values in TablelA.fl 



Table A.2. List of the complementary data. All fluxes have been 
scaled to 80" beam size. Values in brackets were not included in 
the fit. 
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Appendix B: The clumpy PDR model 

To model the far-infrared line emission from S140 we used a 
superposition o f spherical clumps described by the KOSMA- 
T PDR model dRollig et al.L l2006h th at represent an ens emble 
of clumps with a fixed size-spectrum (ICubick et al.L l2008i) . The 
KOSMA-T PDR model simulates a spherical cloud with a radial 
density profile given by 



|http://herscheLesac.esa.int/Docs/SPIRE/htmI/spire-om.html| 



n(r) 



c"* for r < Rjc 

(r/R)-^ for R/c <r<R 
foir>R . 



(B.l) 
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The constant c determines the dynamic range that is covered by 
the power-law density decay. The spectrum of PDR clumps is 
characterized by the clump mass spectrum 



dN 
dM 



aM 



(B.2) 



where the factor a is determined by the total mass of clumps 
within the beam. Mens, and the mass-size relation 



M = CR^ , 



(B.3) 



which implicitly defines the surface density of the individual 
clumps Hs. The factor C is determined by the average ensemble 
density »ens- Gui ded by the clump-decomposition results from 
iHeithausen et alj 01998.) we fix the parameters of the spectrum 
to a = 1.8, 7 = 2.3, 6 = 1.5, and c - 5. The boundaries of the 
clump size distribution M,j,in and M^ax are used as free param- 
eters with the constraint that the maximum clump mass cannot 
exceed half of the mass of the total ensemble. 

Every individual clump is treated as a spherically symmetric 
configuration illuminated by an isotropic external UV radiation 
field, specified in terms of the average interstellar radiation field, 
Xo - 2.7 X 10"^ erg cm"^ s"' in Draine units, and cosmic rays 
producing an average ionization rate, ^cr = 5 x 10"'^ s"' . The in- 
ternal velocity dispersion of molecules within the clumps is fixed 
to 1 km s '. The model computes the stationary chemical and 
temperature structure by solving the coupled detailed balance of 
heating, line and continuum cooling and the ch emical network 
using the UMIST data base of reaction rates dWoodall et al.L 
2007) expanded by separate entries for the '^C chemistry (see 
Rollig et al., 2007, for details). The chemical network currently 
does not include "^O, so that C'^O predictions can only be ob- 
tained by scaling the '^CO values ignoring fractionation between 
the two species. As this ignores the different self-shielding of 
^^CO and C'^O, the model results for C'^O are less rehable than 
for the other lines. 

In the superposition of clumps, the line emission from the 
different clumps is simply added assuming that the velocity dis- 
persion between the clumps is large enough, so that they do not 
shield each other in position-velocity space. This is valid for 
most species, except for the [Oi] emission which is optically 
very thick (r ^ 100), so that the lines are much broader than 
the velocity distribution. Therefore, the model is not able to give 
any reliable estimate for the [Oi] intensities. For the continuum 
extinction of the UV radiation, the situation is different. There, 
mutual shading of the clouds is relevant, leading to the concept 
of different clump ensembles that "see" different UV fields if the 
average ensemble extinction exceeds an Ay of about unity. 

The parameters of the clump mass spectrum imply that most 
of the mass is actually contained in the largest clumps that also 
have the maximum column density or Ay, respectively. The de- 
pendence of Ay and clump size on clump mass is given by 



R^5.3x 10 



\n[cm^^] 



(B.4) 



and 5 from lCubick et al.l (l2008h show that the [CII] emission is 
dominated by the largest clumps while the high-7 CO lines are 
dominated by the smallest clumps. A complex, non-monotonic 
behavior is observed for the lines from atomic oxygen and mid- 
to low-y CO isotopes. 

As the clump spectrum is purely observationally based, it has 
no direct relation to a stability criterion. In fact we find that both 
the most massive clumps are unstable to gravitational collapse 
and the smallest clumps will be dispersed on the timescale of a 
few million years. Therefore, the spectrum can only be consid- 
ered as a snapshot of interstellar turbulence that reflects the den- 
sity structure over a timescale of 10'' - 10^ years. The assump- 
tion of a steady-state chemistry and energy-balance is therefore 
only applicable if all rates are faster. This holds for the dense 
clumps in the S140 model fit with densities above lO'* cm"^, but 
for lower de nsities, an explici tly time dependent modeling would 
be required (IViti et al.Ll2006l) . 



and 



Ay = 1.6018 X 10^^ n[cm-^]R[cm] 



(B.5) 



In terms of the total clump surface or the total solid angle of 
the different clumps we find about equal contributions from each 
logar ithmic mass bin in the ensemble (Eq. 16 in Cubicket al], 
|2008|) . Consequently, we find a non-trivial dependence of the 
intensity in the different HIFI lines on the clump mass. Figs. 3 



